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Introduction 


□ Develop in-house CFD codes to support NASA Armstrong flight research projects 

❖ Hyper X/X-43 project 

❖ Adaptive Compliant Trailing Edge (ACTE) project 

□ Validate the in-house CFD codes 

❖ Compare with other commercially available CFD codes 

❖ Compare with wind tunnel and flight test data 




Summary of Flow Equations 


“NAVIER-STOKES” Eqn (Viscous, turbulent flow) 



DISCRETIZATION TECHNIQUES 


FINITE DIFFERENCE FINITE VOLUME FINITE ELEMENT 
(FD) (FV) (FE) 
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Finite Element Formulation 


□ Navier-Stokes equation: 


dv dfi dg t 
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where 


v = [P pUj pE] T , j = 1,2,3 
fj = [puj ( pUiUj + pStj ) Uj(p + pE)] T , ] = 1,2,3 
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□ Nondimensionalized for numerical calculation, flfy becomes: 

9j = [° (. u i a ij ~ <lj)\ 

And 


^ Re 
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Finite Element Formulation (Cont.) 


□ Taylor’s expansion of the solution v(x, t ) in the time domain: 

A 17 = —At 


dft + dg t 


dxt dx iJ(t) 

□ Applying Galerkin’s spatial idealization 17 = av, the flow equation can be expressed as 

dUi 


M Av = —At 


dxt 


M + K 


v — A t(f 1 +f 2 ) + AtR + At [K a + f a ] 


□ A novel two-step solution procedure is adopted for the flow equation, the inviscid solution being 
augmented with the viscous term and stabilized with artificial dissipation terms. Assuming 

Av = v n+1 - v n 

Then 


At \ At 

1+ 2 C ) M+ 2 K 


Vn+1 


At \ At 

1 ~y c ) m ~t k 


v n AtR 


Or 

Let 


[M+]v n+1 = [M_] v n + AtR 
M+=D++ M'+ 


The viscous stress tensor and heat flux may then be solved as follows: 
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Finite Element Solution Procedure 


□ 

□ 

□ 

□ 

□ 


Step 1: Form 


[D + ]v n+1 = [M.]v n - [M' + ]v n+1 + A tR 


Step 2: Solve v n+1 iteratively 



+ + 
II 
+ 

1 

_1) {[M- ]v„ - [M\ ± + A t(R + R + K a + f a )} 

Step 3: If 

w n+ 1 

* EPS1 

1 ? (0 
t/ n+ 1 

go to Step 2. 

Step 4: If 

17 71 + 1 

* EPS2 

v (i) 

t/ n+ 1 

go to Step 1 . 


Step 5: Repeat Steps 1 to 4 NITER times until desired convergence is achieved, that is until v n+1 ~ 
v n ; EPS1 and EPS2 are suitable convergence criteria factors, specified by the users. 


Note: The iterative process in Step 2 requires a small number of steps, usually 1, and achieves a 
stable, convergent solution. In regions of high pressure gradients, artificial dissipation term is 
applied to prevent oscillations near discontinuities. This is implemented by incorporating pressure- 
switched diffusion coefficients as appropriate. 

This procedure is adopted in the STARS -CFDSOL code that enables effective solution of the 
Naviar-Stokes equation in most flight regimes 
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Numerical Example 


□ 

□ 


Gulfstream III (GUI) Wing finite element discretization 

❖ 15,538 grid 3 1 ,072 elements on wing surface 

♦♦♦ 226,444 grid 1,256,879 million elements for the solution domain 

Flight condition: 

❖ Mach 0.701, 

❖ Angle of attack 3.92 
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Pressure Distributions 




C p plot at span station 145 (368.3 cm spanwise location) 
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Pressure Distributions 



C p plot at span station 230 (584.2 cm spanwise location) 
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Pressure Distributions 



C p plot at span station 385 (1003.5 spanwise location) 
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XXK 


▲ 


C distribution on wing surface 
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C p at station 145 
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C p at station 230 
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C p at station 385 
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CFD Codes Comparison 


CFD Solver 

Flow Equation 

Platform 

No. of CPU 

Total CPU time 

Grid Size 

Note 

STARCCM++ 

RANS, finite 
volume, K-omega 
SST turbulence 

Cluster 

00 

o 

6hrs., 40min (533 
CPU hours) - 3000 
iterations 

7.2M 

polyhedra/prismatic 
for half model without 
T-tail 

number of 
processers is an 
estimate, and the 
time is an 
estimate for that 
number of 
processors 

Fluent 

RANS, Finite 
Volume, K-Omega 
SST Turbulence 

Cluster 

32 

16 hrs., (512 CPU 
hours) 

8.12 M 

tetrahedral/prismatic 
for half model without 
T-tail 

Time is 
approximate 

STARS (MG) 

Euler, finite 
element 

Dell M620 8GB 
Ram, 64 bit 

1 Intel Core i7 
@2.67 GHz 

2.8 hrs. ,(100 
steps, 25 inner 
cycles) 

1.2 M Tetrahedrons for 
wing only 


STARS 

(CFDSOL) 

Full N-S, finite 
element 

Dell M620 8GB 
Ram, 64 bit 

1 Intel Core i7 
@2.67 GHz 

13.8 hrs., (10000 
steps) 

2.8 M Tetrahedrons for 
wing only 


USM3D 

Full N-S, finite 
volume 

Mac 64 bit 

2 CPUs 

16 hrs. 

1.9 M cells for half 
model without T-tail 


FUN3D 

RANS, Finite 
Volume with 
turbulence model 

Cluster 

196 

6 hrs. 

20.4 M nodes for half 
model without T-tail 


TRANAIR++ 

Full potential + 
viscosity 
(boundary layer) 

Linux 

Workstation 

1 CPU 

2hrs., 28min 

1.7M cells for full 
model 
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Concluding Remarks 


□ This paper presents detailed comparison of solutions of the GUI aircraft wing obtained by a number 
of commercially available CFD codes as well as two NASA AFRC in-house codes that use a finite 
element discretization employing unstructured grids. 

□ Finite element formulation of the novel in-house CFD code (CFDSOL) is also presented. 

□ Each of the codes shows reasonable correlation with flight and wing tunnel test data; CFDSOL and 
MG codes appears to be rather close to the two test results. 

□ Use of a single CPU to derive solutions testifies to their cost effectiveness. 
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CFDSOL Flow Chart 
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